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'CJs'- ABSTRACT 

0^ 

In this paper, we argue that in the symbiotic star RR Tel the existence of an accre- 
1-^ ■ tion disk around the hot companion is strongly implied by the characteristic features 

["T , ■ exhibited by the Raman-scattered O VI lines around 6830 A and 7088 A. High degrees 

QQ . of polarization and double-peaked profiles in the Raman-scattered lines and single-peak 

profiles for other emission lines are interpreted as line-of-sight effects, where the H I 
scatterers near the giant see an incident double-peaked profile and an observer with a 
^ ' low inclination sees single-peak profiles. It is predicted that different mass concentra- 

ly-v I tions around the accretion disk formed by a dusty wind may lead to the disparate ratios 

CN ' of the blue peak strength to the red counterpart observed in the 6830 and 7088 features. 

f— ^ ■ We discuss the evolutionary links between symbiotic stars and bipolar protoplanetary 

^N . nebulae and conclude that the Raman scattering processes may play an important role 

in investigation of the physical properties of these objects. 
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1. Introduction 

Symbiotic stars exhibit the thermal components typical of a cool star and a hot star in their 
spectra with additional emission nebulosity. It is usually thought that they form a binary system 
of a giant suffering a heavy mass loss and a white dwarf surrounded by an emission nebula (e.g. 
Kenyon 1986). The double-peaked profiles often observed in the emission lines of many symbiotic 
systems convincingly imply that the emission regions may be characterized by disk-type motions. 
The variability and nova-like outbursts may also be attributed to the existence of an accretion disk. 
Since the giant provides material to the hot star in the form of a stellar wind, it is plausible to 
expect that an accretion disk may be formed around the hot star (Robinson et al. 1994). Recently, 
Mastrodemos & Morris (1998) presented their numerical computations on disk formation in a wide 
binary system of a white dwarf and a giant through a dusty wind. 

An important and distinct aspect of spectroscopy of symbiotic stars is provided by the Raman- 
scattered features around 6830 A and 7088 A, which are identified by Schmid (1989). According 



to him, they are originahy the O VI 1032, 1038 doublet lines that are absorbed and re-emitted 
by atomic hydrogen initially in the ground Is state and finally to de-excite to the 2s state. Basic 
atomic physics of the Raman-scattering processes is discussed by many researchers (Lee & Lee 
1997a, Nussbaumer et al. 1989, Saslow & Mills 1969, Sadeghpour k Dalgarno 1992). 

The Raman-scattered features are characterized by the strong polarization and the Doppler 
enhancement by a factor of ~ 7 attributed to the incoherence of scattering. Therefore, the broad- 
ened profiles and polarization are expected to carry much information about the physical properties 
of the emission regions and scattering geometry. This point is illustrated by the numerical works 
by Harries &: Howarth (1997) and Schmid (1992) and by the spectropolarimetric observations 
performed by Harries & Howarth (1996), and Schmid & Schild (1990). 

However, thus far it appears that neither numerical works nor observational reports successfully 
address the important role of the Raman-scattered features that they provide a unique view of the 
emission region from the neutral scatterers rather than in the line of sight to the observer. One 
good example of this point is found in the spectropolarimetry of the broad double-peaked Ha line 
in the radio galaxy Arp 102B by Corbett et al. (1998). They proposed that the single-peaked 
profile in the polarized flux is mainly associated with the different line of sight to the scatterers 
than the observer's sight line that is prone to give a double-peaked profile in the direct flux. Similar 
but more interesting effects are expected in the case of the Raman-scattered lines, because these 
are purely scattered features with no direct component and furthermore the associated polarized 
fluxes carry additional and complementary information. 

Recently, simultaneous observations of the Raman-scattered features and emission lines in the 
UV range have been performed (Espey et al. 1995, Birriel et al. 1998, Schmid 1998). In this 
regard, this viewpoint about the Raman-scattered lines is expected to be useful for building more 
elaborate models of symbiotic stars. In this Letter, we adopt this point to compute the profiles and 
polarization of the Raman-scattered O VI lines and explain many prominent features seen in the 
symbiotic star RR Tel, and discuss the evolutionary status in relation to bipolar protoplanetary 
nebulae. 



2. Model 

2.1. Accretion Disk in Symbiotic Stars 

In this subsection, we give a brief description of our model adopted for the computation of the 
profiles and polarization of the Raman-scattered lines in the symbiotic star RR Tel. Spectropo- 
larimetric observations of RR Tel have been performed by a number of researchers (e.g. Harries 
& Howarth 1996, Espey et al. 1995). The main features of the spectra of RR Tel include the 
clear double-peaked profiles, high degree of polarization and polarization flip in the red wing in 
the Raman-scattered lines. Because the scatterers responsible for the Raman-scattered lines are 



believed to be located near the giant, the double-peaked profiles imply that the emission regions 
are characterized by a disk-type motion. 

In Fig. 1 is shown a schematic diagram illustrating the accretion disk emission regions and the 
scattering geometry adopted in this Letter. There are two main emission regions around the hot 
star, that is, the red emission region (marked by 'RER') and the blue emission region (marked by 
'BER'), which give the red and blue components in the direction of the binary axis, respectively. 
Here, by the binary axis we mean the line connecting the two stars. According to Mastrodemos & 
Morris (1998) more mass concentration is expected in the RER than in the BER. We discuss this 
effect in the next section. 

The scatterers are assumed to be concentrated mainly in two regions, that is, in region (A) 
near the giant star and in region (B) forming a slowly expanding shell. Solf (1984) investigated the 
bipolar mass outflow of the symbiotic star HM Sge and found that in addition to the outflow of 
velocity ~ 200 km s~^, there are slowly moving features near the equatorial plane. The amount of 
neutral hydrogen in region (B) is not certain at the moment, and in this paper it is assumed to be 
much smaller than that in region (A), but not negligible. It is also assumed that when the incident 
wavevector kj makes an angle less than 45° with the binary axis, then the photon hits region (A) 
and otherwise it is scattered in region (B). 



2.2. Double-Peaked Emission from an Accretion Disk 

We follow the approximation adopted by Home & Marsh (1986) to generate the line profiles 
from the disk inclined at various angles. The emission lines are assumed to be optically thick, and 
the disk obeys the Keplerian rotation. We also assume that the disk is geometrically thin with the 
disk height ratio H/R = 1/50. Since we do not know enough about the accretion disk in symbiotic 
stars and even less about the emission regions in the disk, we just take the line source function 
to be a simple power-law Sl oc r^"*^'^, which has been used by Home & Marsh (1986) to match 
the line shape of Ha in Z Cha. Fig. 2 shows the line profiles for inclination angles i = 10°, 30°, 
60°, 90°. The maximum velocity of the emission region is chosen to be 50 km s~^ and the thermal 
broadening has not been applied. We note that the low-inclination profile is narrower than the 
high-inclination counterpart. 

In this paper, we do not attempt to reproduce the observational results by fitting accurately, 
and instead we want to point out the physical origins of the main characteristic features shown in 
the observational data. Therefore, for the computation of the generic profiles and polarization, we 
use two Gaussians to approximate the profiles in Fig. 2, and adopt the peak strength ratio as a 
new free parameter. It is referred to Lee & Lee (1997b) for the detailed description of the Monte 
Carlo computation of the scattered flux profile and polarization. (See also Harries & Howarth 1997, 
Schmid 1992) 



3. Result and Discussion 

3.1. Profile and Polarization 

In Fig. 3 is shown the main result. We set the strength ratio of the blue part to that of the 
red part to be 2/3 and discuss this point in the next subsection. The top panel shows the flux and 
polarization of the component scattered near the giant. In Panel (b) the same quantities are shown 
for the radiation scattered in the spherical shell receding with velocity Vsheii = 30 km s~^ and total 
scattering optical depth tt = 0.3. In the bottom panel is shown the sum of the preceding two 
components. Here, the polarization direction is represented by the sign, where the positive sign 
represents the polarization in the direction perpendicular to the binary axis whereas the negative 
sign the parallel direction. 

The double-peaked profile in the top panel is obtained because the scatterers near the giant 
see a double-peaked incident source. This component is strongly polarized in the direction perpen- 
dicular to the binary axis. Because the relative velocity of the scatterers near the giant is ignored, 
the scattered flux profile is almost the same as the incident profile in the binary axis direction and 
the breadth of the feature represents mainly the kinematics of the source part. 

In Panel (b), we obtain a broad profile with a single peak. The polarization is in the direction 
parallel to the binary axis and is weaker because region (B) subtends a larger solid angle than 
region (A) does. Due to the relative motion of the shell, the location of the peak is redshifted from 
the center. The single-peaked profile is obtained because the averaged profile incident to the shell 
is smooth with a single peak near the center. 

In combining the two components we may see interesting fine points that may be revealed 
in the spectropolarimetry of symbiotic stars. Firstly, the strength of the scattered component is 
determined mainly by the scattering optical depth and the solid angle of the scattering region. Our 
choice is such that the shell has a small scattering optical depth so that the synthetic scattered flux 
is dominantly determined by the near axis-scattered component. In RR Tel, most emission lines 
are single-peaked and the Raman lines show the double-peaked profile, which strongly imply the 
emission regions in disk-type motion with low inclination. 

Secondly, the shell-scattered component is red-shifted and single-peaked. Therefore, this com- 
ponent adds more flux to the red part and increases slightly the ratio of the red peak strength 
to that of the blue in the synthetic flux. However, since the shell-scattered component is polar- 
ized in the parallel direction, in the overlapping region with the flux scattered in region (A) the 
total polarized flux is reduced as a result of cancellation of polarization but remains still in the 
perpendicular direction. Because the shell-scattered component extends more redward than the 
component scattered in region (A), there remains a parallelly polarized flux in the reddest part 
where the (A)-scattered component does not contribute. The polarization flip in the red wing part 
is an important feature shown in many spectropolarimetric observations of symbiotic stars. 



Finally, as Harries & Howarth (1996) pointed out, the overall profiles observed for many 
symbiotic systems are broader than the terminal speed of the stellar wind associated with K or M 
type giants and the commonest profile type is triply-peaked where the reddest component is often 
polarized oppositely compared with the remainder part. It turns out that the overall breadth of the 
profile is determined by the kinematics of the source part and the relative motion of the scatterers. 
Furthermore, by increasing the speed of the receding shell, we can also generate triply-peaked 
profiles with the polarization flipped in the reddest component. 



3.2. O VI Doublet Strength Ratio 

In this subsection we discuss a possible prediction of the symbiotic star model consisting of 
an accretion disk with a bipolar wind. According to Mastrodemos &: Morris (1998), it is expected 
that the RER has a higher optical depth than the BER does. For the O VI 1034 doublet, the 
1032 photons have twice larger optical depth than the 1038 counterparts. In an optically thin 
medium, the oscillator strength ratio becomes that of the emergent line strength ratio. However, 
in an optically thick medium the emergent line strength of the 1032 A component becomes similar 
to that of the 1038 A component. However, complications may occur in a non-stationary medium, 
where the sensitive dependence of the resonance scattering cross section on the velocity field may 
easily alter the escape probability so that it becomes eventually proportional to the velocity gradient 
(e.g. Sobolev 1963, Lee & Blandford 1997, Michalitsianos et al. 1988). Further complications can 
be expected if the medium is dusty, in which case resonantly scattered O VI photons are subject 
to destructions by the dust particles (e.g. Neufeld 1991). 

It is an interesting possibility that if the BER is optically thin and the RER is thick, then the 
Raman-scattered line around 6830 A will show a larger ratio of the blue peak strength to the red 
peak strength than the 7088 A feature, as is depicted in Fig. 1. It appears that in RR Tel the 
7088 A feature has a weaker red part compared that of the 6830 A feature (Harries & Howarth 
1997). However, since the 7088 A feature is much weaker, higher resolution spectroscopy with good 
signal to noise ratio must be invoked to exclude other possibilities such as selective interstellar 
absorptions. 

Hayes & Nussbaumer (1986) proposed that the electron density ne ~ 3 x 10^ cm~^ and the 
size of the emission line region R ^ 10^^ cm in RR Tel. A simple computation of the line center 
optical depth of the O VI 1034 doublet in the emission region gives 

Te ~ 7.2 X 10^T7'/'[/,^ov//10-']K/(106 cm-3)][i?/(100 AU)], (3.1) 

where T4 is the temperature in units of 10^ K, Aqvi is the O VI fraction in number and fi is the 
oscillator strength of the O VI resonance transition (Rybicki &: Lightman 1979). Apparently, the 
emission region in RR Tel is very optically thick and the O VI doublet strength ratio would be 
nearly 1:1 in the entire region, even if the BER is an order of magnitude rarer in mass concentration 



than the RER. Therefore, both high resolution spectroscopy and further theoretical work on the 
radiative transfer based on a more refined model will shed light on this point. 



3.3. Evolutionary Status of Symbiotic Stars 

Asymmetric morphologies in the emission nebulae have been known in many symbiotic systems 
including V 1016 Cyg (Schild & Schmid 1996), HM Sge (Solf 1983, Eyres et al. 1995). Solf(1984) 
emphasized that a large fraction of symbiotic stars exhibit bipolarity and that there is a remarkable 
morphological resemblance with postnova shells and protoplanetary nebulae. The spectroscopic 
similarity between symbiotic stars and bipolar protoplanetary nebulae has led many researchers 
to propose evolutionary links between them (e.g. Iben & Tutukov 1996, Corradi 1995, Corradi & 
Schwarz 1995). A supporting argument to this effect is that the binarity may play an important role 
in forming an aspherical morphology of the nebula, typically characterized by bipolarity (Morris 
1987, Soker 1998). However, the binarity of symbiotic stars is well-established whereas that of the 
bipolar protoplanetary nebulae still remains controversial. 

The accretion disk formation in a young planetary nebula has been discussed by a number 
of researchers (e.g. Morris 1987, Soker & Livio 1994). Recently, Mastrodemos &: Morris (1998) 
presented an interesting numerical computation on the dusty wind accretion in a detached binary 
of a mass-losing AGB star and a hot star that may be responsible for the bipolar morphology in 
protoplanetary nebulae. They found that a permanent and stable accretion disk is formed around 
the hot companion with an efficient cooling associated with dust in the wind (see also Theuns & 
Jorissen 1993). In particular, they concluded tha the limiting binary separation for disk formation 
should be greater than 20 AU for their M4 model. Therefore, it is interesting that the evolutionary 
links between symbiotic stars and bipolar protoplanetary nebulae imply that the hypothesis of 
accretion disk emission in RR Tel is also rendered a strong case. Furthermore, as shown in Fig. 1 
the bipolar wind along the disk axis may provide natural scattering sites for the Raman-scattered 
flux that constitutes the red wing part with polarization flip. 

The Raman-scattered lines are characterized by the high polarization and broadened profiles, 
which enable one to put strong constraints on the scattering geometry. Therefore, the main spectro- 
scopic and polarimetric features imply very plausibly that the symbiotic star RR Tel may possess 
an accretion disk with a bipolar outflow. So called 'type 2' symbiotic stars seem to show simi- 
lar behaviors in the Raman scattered lines and exhibit bipolar type morphologies. Furthermore, 
about a half of symbiotic stars exhibit the Raman scattered lines. Fine-tuned conditions such as 
a restricted range of the initial mass ratios of the two constituent stars and the orbital parameters 
affecting the mass transfer rate may be needed for the co-existence of a large amount of neutral 
hydrogen and highly ionized nebulae which characterize the symbiotic phenomenon. 

An interesting example is provided by Pequignot (1997), who performed high resolution spec- 
troscopy on a young planetary nebula NGC 7027 and found a Raman scattered He II line blueward of 



H/3. Gurzadyan (1996) discusses the relation of NGC 7027 to symbiotic stars, noting that NGC 7027 
shows high ionization hnes along with a strong IR component. So far the Raman-scattering by H I 
is found to operate only in symbiotic stars except NGC 7027, and theoretical possibilities have been 
discussed in other astrophysical objects such as active galactic nuclei (Nussbaumer et al. 1989, Lee 
&: Yun 1998). Therefore, it will remain an interesting possibility that similar processes are expected 
to operate in bipolar protoplanetary nebulae, in which case the Raman scattering can be regarded 
as an important tool that reveals the evolutionary links between symbiotic systems and bipolar 
protoplanetary nebulae. 
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Fig. 1. — A schematic diagram illustrating the accretion disk emission region and the scattering 
geometry for the symbiotic star RR Tel adopted in this paper. Conceptually there are two emission 
regions around the hot star denoted by 'RER' and 'BER', which provides red-shifted line photons 
and blue-shifted ones to the direction of the giant respectively. The scattering region is largely 
divided into two regions, i.e., region (A) near the giant that is very thick and region (B) that is 
a large spherical shell expanding slowly with speed Vsheii = 30 km s~^ and total scattering optical 
depth tt = 0.3. The scatterers near the giant see a double peak incident source whereas the source 
appears to have a single-peaked profile. It is assumed that when the incident wavevector kj makes 
an angle less than 45° with the binary axis, then the photon hits region (A) and otherwise it enters 
region (B). 
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Fig. 2. — The emission line profiles from an optically thick accretion disk for inclination angles 
i = 10° (solid line), 30° (dotted line), 60° (dashed line), and 90° (long-dashed line) obtained by using 
the approximation adopted by Home & Marsh (1986). Note that i = 90° corresponds to the edge- 
on disk case. The maximum velocity of the emission region is chosen to be 50 km s^^ and the 
thermal broadening has not been applied. 
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Fig. 3. — Profiles (dotted lines) and polarization (solid lines with la error bars) of the Raman- 
scattered O VI lines from accretion disk type emission regions shown in Fig. 1. In Panel (a), the 
flux and polarization of the component scattered in region (A) of Fig. 1 are shown. In Panel (b), 
the same quantities as in Panel a for the component scattered in region (B) are shown. In Panel 
(c) are shown the synthetic profile and polarization of the quantities shown in Panels (a) and (b). 



